Reconstruction of nuclear charged fragment 
trajectories from a large gap sweeper magnet 



N. Frank a ' b ' c '*, A. Schiller a , D. Bazin a , W.A. Peters a ' b , and 

M. Thoennessen a ' b 

^National Superconducting Cyclotron Laboratory, Michigan State University, East 

Lansing, MI 48824, USA 

h Dept. of Physics & Astronomy, Michigan State University, East Lansing, MI 

48824, USA 

c Dept. of Physics, Concordia College, Moorhead, MN 56562, USA 



Abstract 

A new method to reconstruct charged fragment four-momentum vectors from mea- 
sured trajectories behind an open, large gap, magnetic dispersion element (a sweeper 
magnet) has been developed. In addition to the position and angle behind the mag- 
net it includes the position measurement in the dispersive direction at the target. 
The method improves the energy and angle resolution of the reconstruction signif- 
icantly for experiments with fast rare isotopes, where the beam size at the target 
position is large. 
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1 Introduction 

The frontier of nuclear physics is moving to more and more neutron rich 
systems extending out to the limit of existence (dripline). Only a few if any 
bound excited states exist in nuclei close to the neutron dripline and beyond, 
none at all. The unbound state(s) of nuclei near the neutron dripline are 
energetically above one- or even two-neutron separation energies. Thus these 
nuclei break-up by neutron emission, which makes gamma spectroscopy an 
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ineffective tool for the study of these nuclei. Neutron unbound states have to 
be reconstructed by neutron decay spectroscopy [1-4]. 

Often, the neutron-rich nuclei of interest are produced via one or more nucleon 
stripping reactions from rare-isotope beams produced in fast fragmentation 
reactions. The lifetime of these systems is extremely short (~10~ 21 s) so they 
decay within the target immediately after production. Nuclei produced by the 
fast fragmentation technique have high energy (~100 MeV/A) and intensity 
(10 4 /s) and, beams of such nuclei in general have a large spread in momentum, 
size, and angle. Due to the high energy of the beams, the charged fragments 
as well as the neutrons are extremely forward focused and have to be detected 
around zero degrees. The decay energy from neutron unbound states can be 
reconstructed from a measurement of the four-momenta of the decay products 
(charged fragment and neutron). 

One method of detecting both the charged fragments and neutrons in coinci- 
dence is to separate the charged fragments from the neutrons using a dipole 
magnet. In order to maximize the angular acceptance of the neutrons at zero 
degrees the dipole (or sweeper) magnet has to have a large gap, an opening 
at zero degrees ("C"-magnet), and has to be located as close to the target as 
possible. These stringent requirements do not allow the use of a full magnetic 
spectrometer which is typically used to determine the energy and angles of 
emitted particles in nuclear structure experiments [5-11]. For these spectrom- 
eters ion-optical codes have been developed to extract the four-momentum 
of the fragments at the target from the measured quantities following the 
spectrometer. One example for such a code which is extensively used is COSY 
infinity [12-14]. However, as we will discuss below, the standard application 
of COSY INFINITY to the specific problem of track reconstruction for the case 
of a simple, open, large-gap magnet and a large beam-spot size, which are 
typical for neutron-decay experiments, does not yield satisfactory results in 
terms of resolution. 



2 Method 

2.1 Fully Inverse Ion-Optical Matrices 

A forward ion-optical matrix relates coordinates (x, 9 X , y, 9 y , 5)^ at the 
position of the reaction target in front of the sweeper magnet to coordinates 
(x, 9 X , y, 9y, AL)( D ) at the position of the charged particle detector. In this 
work, x is the dispersive and y is the non-dispersive direction of the sweeper 
magnet. The angular coordinates correspond to ratios of transverse to total 
momentum in their respective directions, e.g., 9 X = p x /po- L is the track 
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length of a charged particle along the central trajectory from the reaction 
target position to the detector position. Thus AL is the difference in the track 
length compared to the length of the central trajectory. The central trajectory 
is defined by (x, 9 X , y, 9 y , 5)^ = and (x, 9 X , y, 9 V Y D ^ = 0. The relative 
energy deviation 5 is given by 



E — Eq 
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(1) 



where E is the energy of the incoming particle and Eq is the reference energy 
of a particle following the central trajectory. This energy is determined by the 
measured magnetic field of the sweeper magnet. 

In first order, coordinates belonging to the dispersive direction and the non- 
dispersive direction do not depend on each other because of the magnetic 
field's symmetry in the non- dispersive direction. Thus, a typical first-order 
forward ion-optical matrix will look like the following 
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where there are two non-zero 2x2 sub-matrices M' . The sub-matrices mix 
position and angle in either the dispersive or non-dispersive direction, respec- 
tively. For mixing position and angle in the same direction, det M' = 1 is 
required in order to preserve the phase-space volume P~| In order to adequately 
describe the sweeper magnet's field, the higher-order dependence of output 
coordinates on combinations of input coordinates is necessary. For this pur- 
pose, the input vector is extended to all possible higher-order combinations 
and the ion-optical matrix is modified accordingly. In first order, there are 
only five linear input coordinates, hence the ion-optical matrix is a 5 x 5 ma- 
trix. In second order, there are 15 new quadratic input coordinates such as 
((x^) 2 , x < - T - ) 6 1 ^); the ion-optical matrix has to be extended by a 5 x 15 ma- 
trix. In third order, there are 35 new third-order input coordinates and the 
ion-optical matrix has to be extended by a 5 x 35 matrix and so on. The 
ion-optical matrices produced and used in this work are up to third order. 



1 With the given beam energies and intensities, non-conservative effects from space 
charges and photon production are estimated to be negligible. 
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In typical spectrographic applications, the position and angle of particles at 
the charged particle detector position are used to determine the energy, angle, 
and position of those particles at the reaction target. For this reason, the 
forward ion-optical matrix has to be inverted. However, looking at Eq. (2), 
one notices that a full inversion is not useful because typically only four out 
of the five coordinates of the particle are known behind the spectrograph 
(the track length is not a measured quantity). This fundamental problem 
has been solved in two different ways: (i) if focusing magnetic elements are 
available, the spectrograph can be run in dispersion-matched mode, in which 
essentially the matrix elements M relating the input coordinate x^ to any 
of the output coordinates are negligibly small [14], or (ii) the spectrograph 
can be run in focused mode, in which the incoming beam is focused in the 
dispersive direction at the target, such that one can assume = 0. Either 
way, Eq. (2) can effectively be reduced to 
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where the column concerning the target coordinate x' ' and the row concern- 
ing the track length variation AL have been eliminated. After inverting this 
matrix, the coordinates (9 X , y, 9 y , 5)^ can be reconstructed from position and 
angle measurements by the charged particle detectors. The program COSY 
provides this kind of inverse ion-optical matrix. In this work, these matrices 
will be called fully inverse ion-optical matrices. 

If it is experimentally not feasible to run in dispersion-matched or focused 
mode, the above assumption of x^ = is not a valid assumption. The use 
of the fully inverse ion-optical matrix in the standard COSY method will not 
result in sufficient resolution of the reconstructed four-momentum. In these 
cases it is necessary to measure the beam position x^ at the target and 
include this information in the reconstruction procedure. 

2.2 Partial Inverse Ion-Optical Matrices 

The reaction target coordinates x^ and y^ can be measured with tracking 
detectors located in front of the target. With these two additional pieces of 
information, the problem of reconstruction becomes overdetermined. In the 
following, the measured target coordinate x^ is used as an additional input 
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and entered into the reconstruction. 



The starting point for our discussion is the first-order ion-optical matrix which 
relates only the coordinates in the dispersive direction 
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Here, the target coordinate x^ is an input to the matrix multiplication on 
the right-hand side, and the detector coordinates x^ and 9 X D ^ are outputs on 
the left-hand side. All three are measured quantities. The unknown quantities 
are the target coordinates 9 X T ^ and 6 on the right-hand side and the deviation 
in track length from the length of the central trajectory AL on the left-hand 
side. It is now straightforward to exchange one coordinate on the left-hand side 
with one coordinate on the right-hand side by using the rules shown in Table 
1. The choice of the coordinates which one would like to exchange determines 
the pivot, i.e., the matrix element which relates the two coordinates. 

The ion-optical matrix 
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is generated by exchanging the reaction target coordinate 5 with the detector 
coordinate x^ D \ In a second step, one can now exchange the reaction target 
coordinate 9 X T ^ and the detector coordinate 9 X D ^ using the same rules, such 
that the three known quantities are on the right-hand side and the three 
unknown quantities are on the left-hand side. 

Ion-optical matrices, where coordinates have been exchanged in this fashion 
will be called partial inverse ion-optical matrices in this work. It is straight- 
forward to generalize this method to the 5x5 matrices of first order. At this 
point, we would like to make three comments: (i) exchanges of two coordinates 
are only possible if the pivot does not equal zero, i.e., it is only possible to 
exchange coordinates in the dispersive plane and in the non-dispersive plane 
among themselves, (ii) numerically, the method works better the larger a pivot 
element is chosen. Ideally, the pivot element that is chosen should be greater 
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than one. (iii) when exchanging all five coordinates along the diagonal, the 
resulting ion-optical matrix is the mathematical inverse of the original ion- 
optical matrix, i.e., in first order, the result of the inversion method is math- 
ematically exact. 

The pertinent rules for exchanging coordinates in second or higher order can 
be derived from the following considerations. Assume the input, i.e., target co- 
ordinates of a forward ion-optical matrix are (xj) and the output, i.e., detector 
coordinates are (yj) with i,j = 1...5, then the coefficients of the ion-optical 
matrix are ag^^^^, such that 



yj=1...5 — 2^ a n 1 ,n 2 ,m,n 4 ,n 5 X l x 2 x 3 X 4 X 5 \ U ) 

l<ni+n2+ri3+n4+nr,<N 

where N is the order of the ion-optical matrix and n(j =i ... 5 ) > 0. Without loss 
of generality, consider only the exchange of the target coordinate x 5 with the 
detector coordinate y 5 . Then, the matrix elements n5 defined by 
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which relate the new input coordinates (xj, 1/5) with i = 1...4 to the new 
output coordinates (yj, x 5 ) with j = 1...4 are of interest. In the first step, 
Eq. (6) for j = 5 is entered into Eq. (8) and, order by order, a system of M 
linear equations is constructed from which the coefficients are determined 
by simple matrix inversion. In the second step, Eq. (8) is entered into Eq. (6) 
for j = 1...4 and by a comparison with Eq. (7) in a similar fashion as in the 
first step, the coefficients M 1 --- 4 ) can be deduced. 

Table 2 shows the dimension M of the system of linear equations which has 
to be solved when exchanging one coordinate. The dimension M is given as 
function of the order of the ion-optical matrix (assuming five input and output 
variables in the ion-optical matrix). 

A program has been developed using the algorithm for exchanging coordi- 
nates in ion-optical matrices of five input and output coordinates up to third 
order as described above. Using this program, forward ion-optical matrices 
were transformed into inverse ion-optical matrices. Since the position is 
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available as an input, and 9P are exchanged with 9P and 5^ T \ In ad- 
dition, and 9[ D) are exchanged with y^ and 9^\ After the exchanges, 
the partial inverse ion-optical matrix is obtained as 
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The final result can be checked against the fully inverted ion-optical matrix 
provided by COSY, because reducing the matrix by the fifth row and the fifth 
column should produce the same result as COSY coefficient by coefficient. It is 
important to note that the first exchange of two coordinates in the dispersive 
as well as the non-dispersive plane should be done for a pair of coordinates 
which shares a large (ideally greater than one) first-order matrix element (the 
pivot). 



3 Theoretical Comparisons 



In the partial-inversion process, forward ion-optical matrices with terms up 
to third order are transformed into inverse ion-optical matrices. The inverse 
ion-optical matrices also have terms up to third order. However, a completely 
inverse ion-optical matrix would have fourth and higher order terms. Before 
comparing the standard COSY fully inverted ion-optical matrix to the one 
obtained by the partial-inverse method, the impact of not including these 
higher order terms is quantified. 

Figure 1 shows comparisons of input (thick lines) to output (thin lines) con- 
tours of the relative energy versus target angle using the partial-inverse method 
(top) and the COSY method (bottom) for three different x^ positions. The 
output contours are produced by propagating the input contour through a 
forward transformation and then an inverse transformation using either the 
partial-inverse or the COSY fully inverted ion-optical matrix. The input con- 
tours correspond to a spread in (±3%) and angle 9P (±5 mrad) consistent 
with the characteristics of rare-isotope beams produced in fast fragmentation 
reactions. The top row shows that the differences due to the lack of fourth 
and higher order terms for the partial-inverse method are very small. How- 
ever, with the exception of the contour pair at cm in x^, the comparisons 
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Fig. 1. Contours of relative energy versus target angle. The input contours (thick 
lines) are compared to output contours (thin lines) for the partial-inverse method 
(top) and the COSY method (bottom). The three columns are for positions (in 
cm) of —1 (left), (center), and 1 (right). 

show large deviations for the COSY method. It should be pointed out that 
the difference between input to output contours increases with the size of the 
contour regardless of the method chosen. In general, any transformation not 
performed up to all orders will always show some deviation. 



4 Comparison with Data 



The two different reconstruction methods were compared to actual data from 
an experiment performed at the National Superconducting Cyclotron Labo- 
ratory (NSCL) at Michigan State University. A primary beam of 140 MeV/u 
40 Ar impinged on a 893 mg/cm 2 Be production target. The desired secondary 
beam of 26 Ne at 86 MeV/u was purified and delivered to the experimental 
setup using the A1900 fragment separator [15]. The secondary beam inter- 
acted with the reaction target at the experimental setup resulting in charged 
fragments and neutrons. The charged particles were bent out of the way us- 
ing a large gap sweeper magnet built at the National High Magnetic Field 
Laboratory (NHMFL) at Florida State University (FSU) [16]. The neutrons 
continued forward at zero degrees and were detected with the Modular Neu- 
tron Array (MoNA) [17, 18]. The experiment was a 2p-n stripping reaction 
from a 26 Ne secondary beam to neutron unbound states in 23 as described 
in [19,20]. 

The beam position at the target was determined by propagating a measured 
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Table 1 

The rules of coordinate exchange for a first order matrix. 
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Table 2 

Dimension M of the system of linear equations which has to be solved when ex- 
changing one coordinate. M is a function of N, which is the order of the ion-optical 
matrix. Five input and output variables are assumed. 
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Fig. 2. Dispersive position distribution x^ of beam particles at the reaction target. 
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Fig. 3. Reconstructed versus measured target angle 9 X (top) and reconstructed 
versus measured energy (bottom) distributions for the COSY method (left) and the 
partial-inverse method (right). The diagonal lines represent a one-to-one correspon- 
dence. 



trajectory through a quadrupole triplet ion-optical matrix. The measured tra- 
jectory was determined by two position detectors before the quadrupole triplet. 
The beam position and size in the dispersive direction is shown in Fig. 2. The 
beam was offset from the center by approximately 3 mm and had a width of 
20 mm. 

The relative energy 5 and the target angle 6^ provide the best comparison 
between both reconstruction methods because they have matrix elements in 
first order that depend on x^ T \ Data were taken without a target so that 
the reconstructed energy could be compared directly with the energy of the 
incoming 26 Ne beam which was determined from a time-of-flight measurement 
between two scintillators located in the beamline separated by 35.7 m. 
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Figure 3 shows a comparison of the reconstructed versus measured target angle 
#( T ) (top) and reconstructed versus measured energy (bottom) distributions 
for the COSY method (left) and the partial-inverse method. The distributions 
from the COSY methods show significantly larger deviations from the ideal one- 
to-one correspondence indicated by the solid diagonal line than the partial- 
inverse method. 

A more quantitative comparison is shown in Figs. 4 and 5 where the differ- 
ences between the measured and the reconstructed target angle and energy, 
respectively, are plotted. The centroid of 9^ distribution is not reproduced 
by the fully inverse COSY matrix and differs by ~6 mrad. The improvement 
for the partial-inverse method is obvious as the centroid of the distribution 
agrees with the measured data. In addition, the result of the partial-inverse 
method yields a narrower width for the difference distribution of 9%F\ 

Even more important is the improvement for the reconstruction of the energy. 
Figure 5 shows the deviation of the reconstructed energy from the measured 
energy in %. On average, a ~2% deviation in beam energy is observed for 
the COSY method, while the partial-inverse method reconstructs the energy 
correctly. Also, the FWHM of the distribution using the partial-inverse method 
is only ~1.5% as compared to ~4% using the COSY method. 
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Fig. 5. Difference between the measured and reconstructed beam energy at the 
reaction target position for the partial-inverse (solid) and COSY (dashed) methods. 

5 Conclusion 

The partial-inverse method significantly improves the accuracy of the recon- 
struction of four-momenta for experiments using secondary beams which have 
typically large beam spots. The assumption of an infinitely small beam spot 
size used in the fully inverse COSY method is insufficient to extract spec- 
troscopy information from the data. The partial-inverse method utilizes the 
information of the finite beam spot size for the correct reconstruction. Thus, 
in order to extract the optimum resolution by using this method it is necessary 
to measure the beam position with appropriate tracking detectors. 
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